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Abstract. We performed a reverse process of a usual optical data analysis of boson- 
exchange superconductors. We calculated the optical self-energy from two (MMP 
and MMP+peak) input model electron-boson spectral density functions using Allen’s 
formula for one normal and two (s- and d- wave) superconducting cases. We obtained 
the optical constants including the optical conductivity and the dynamic dielectric 
function from the optical self-energy using an extended Drude model, and finally 
calculated reflectance spectrum. Furthermore to investigate impurity effects on optical 
quantities we added various levels of impurities (from the clean to the dirty limit) in 
the optical self-energy and performed the same reverse process to obtain the optical 
conductivity, the dielectric function, and reflectance. From these optical constants 
obtained from the reverse process we extracted the impurity dependent superfluid 
densities for two superconducting cases using two independent methods (the Ferrel- 
Glover-Tinkham sum rule and the extrapolation to zero frequency of —ei (w)w 2 ); we 
found that a certain level of impurities is necessary to get a good agreement on 
results obtained by the two methods. We observed that impurities give similar effects 
on various optical constants of s- and d -wave superconductors; the more impurities 
give the more distinct gap feature and the less superfluid density. However, the s- 
wave superconductor gives the superconducting gap feature more clearly than the rf- 
wave superconductor because in the d -wave superconductors the optical quantities are 
averaged over the anisotropic Fermi surface. Our results also supply helpful information 
to see how characteristic features of the electron-boson spectral function and the 
s- and d -wave superconducting gaps appear in various optical constants including 
raw reflectance spectrum. Our study also may help to understand the usual optical 
analysis process thoroughly. Further systematic study of experimental data collected 
at various conditions using the optical analysis process will help to reveal the origin of 
the mediated boson in the boson-exchange superconductors. 
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1. Introduction 

Strongly correlated electron systems including high temperature superconductors 
have been very intriguing research objects and been studied intensively by various 
experimental techniques and theoretical methods [1, 2], Optical and/or infrared 
spectroscopic techniques have been playing an important role to delve into the 
origin of the correlations among electrons in the correlated systems [3]. In general 
measured reflectance spectra of correlated electron systems including boson-exchange 
superconductors [4] can be analyzed, as follows: First a Kramers-Kronig analysis[5] is 
performed to obtain the optical constants including the optical conductivity. Then 
an extended Drude model is applied to get the optical self-energy [6] which can carry 
information of the correlations and is in one-to-one correspondence with the well-known 
quasiparticle self-energy [7, 8, 9]. An extended Allen’s formulas[10, 11, 12,13] are used for 
getting more fundamental quantity, the electron-boson spectral density function, from 
the optical self-energy. The Allen’s formulas relate linearly the optical self-energy to the 
electron-boson spectral density function. The electron-boson spectral density function 
can be obtained usually by solving numerically [14, 15, 16] the extended Allen’s formula 
which is an integral equation. Further systematic optical studies provide temperature 
and doping dependent properties of the electron-boson spectral density function of high 
temperatures superconductors [17, 18]. 

It has been known that the electron-boson spectral functions of cuprates obtained 
from various spectroscopic experimental [angle resolved photoemission (ARPES)[19], 
(scanning) tunneling (STS) [20, 21], inelastic neutron scattering (INS), Raman[22], and 
infrared (IR)[17, 18, 23]] techniques consist of two components: one broad background 
and a sharp mode. While the background shows small doping dependence the sharp 
mode shows strong doping and temperature dependencies [2], The sharp mode appears 
and grows with decreasing temperature by consuming the spectral weight of the 
high temperature broad background[2, 6, 17, 18, 20, 21, 22, 24, 25, 26, 27, 28, 29]. 
The sharp mode is known to be associated with the well-known magnetic resonance 
mode [30]. As reducing the sample temperature the sharp mode starts to appear 
above the superconducting transition temperature (T c ) in underdoped cuprates and 
at T c in the optimally and overdoped cuprates[2, 6, 24, 25, 26]. Many researchers 
believe that the high temperature superconductors including cuprates are boson- 
exchange superconductors as the conventional phonon-mediated superconductors and 
the electron-boson spectral functions can carry information of the bonding force of 
Cooper pairs in the high temperature superconductors. So the electron-boson spectral 
density function observed in cuprates has been known as the superconducting pairing 
glue function[6, 20, 23, 31, 16, 32, 33, 34], Figuring out the microscopic origin of the 
boson which involves in the bonding force is one of the most important problems in 
contemporary condensed matter physics. We believe that the electron-boson density 
functions obtained by various experimental techniques will play a very important role 
to solve the problem. 
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In this paper we performed a reverse process of the usual optical analysis process 
described above to understand better the usual process and investigate impurity effects 
on optical quantities. We started from an impurity-free pure (or ideal) system, and then 
introduced elastic impurities into the system, and finally examined the impurity effects 
on the optical constants including reflectance. To achieve this result we introduced 
two input model electron-boson spectral density functions, which are modeled with 
general shapes of the electron-boson spectral density functions observed by various 
spectroscopic experiments [2], and then calculated the imaginary part of the optical self¬ 
energy (or the optical scattering rate) from the input spectral density functions using 
the Allen’s formulas[10, 15]. One of the two input model spectral density functions 
is the phenomenological antiferromagnetic spin fluctuation [35], which is proposed by 
Millis, Monien, and Pines (MMP) and is also known as an MMP mode. The other 
is a combination of the MMP mode and a sharp Gaussian mode which is modeled 
of an observed electron-boson spectral function in cuprate system at low temperature 
below the superconducting transition temperature [2, 17]. We also considered for three 
different cases: one normal and two different (s- and d- wave) types of superconducting 
states to investigate superconducting gap symmetry effects. Furthermore we used a 
Kramers-Kronig relation to obtain the real part of the optical self-energy from the 
imaginary part (or the optical scattering rate). We calculated the optical conductivity 
using the extended Drude formalism[6, 36]. We also calculated the dynamic dielectric 
function using the relations between optical constants and then finally obtained the 
normal incident reflectance spectrum using the Fresnel equations. 

Impurities are not avoidable experimentally, especially for doped materials with 
substitution. The impurity effect on physical quantities is an important physical issue 
especially for unconventional superconductors including d- wave superconductors [37, 38]. 
Most experimental results may be contaminated by the impurities in the material 
systems or sometimes impurities are introduced experimentally to study their effects 
on the physical quatities[39]. There may be physical issues, which ones may need to 
be aware of impurity effects, for analysis of optical spectra. For example violation of 
the Ferrel-Glover-Tinkham (FGT) sum rule[40, 41] in certain superconductors can be 
interpreted as the lowering kinetic energy which comes from contributions of both the 
condensation energy and the potential energy increase due to larger Coulomb repulsion 
in the Cooper pairs[42, 43, 44], So we added several different impurity scattering rates 
(from the clean to the dirty limit) to the impurity-free optical scattering rate. We 
performed the same reverse procedure to obtain the optical conductivities, the dynamic 
dielectric constants and reflectance spectra from the optical self-energies which contain 
impurities. We discuss impurity effects on other optical constants including reflectance. 
Our results show that measured reflectance spectrum and other optical constants clearly 
show important characteristic features of the input electron-boson spectral function 
and the superconducting gap feature. This may be quite useful because one can see 
important features even in raw reflectance spectra without going through a full optical 
analysis process. And then we investigate impurity effects on the physical quantities 
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including the partial sum, spectral weight redistribution by the Holstein process, and the 
superfluid density which can be obtained from the optical conductivity and the dielectric 
function. Our results show that impurities may affect on the FGT sum rule; one may 
need enough level of impurities in the material systems to hold the FGT sum rule. 

We observed that the absolute superfluid density depends on the impurity level; the 
more impurities give the less superfluid density. We also compared our results obtained 
by the reverse process with experimental data of optimally doped Bi 2 Sr 2 CaCu 2 0 8+( 5 
and observed good agreements for all optical constants considered including reflectance. 
These good agreements confirm that our reverse process is reliable. We could understand 
thoroughly the usual optical process through this study. We believe that further study 
in this direction may help to reveal the microscopic origin of the exchange boson in the 
high temperature superconductors as in the conventional superconductors by collecting 
more experimental data and applying the optical analysis process to them to obtain 
various physical properties of the exchange boson. 

This paper is structured with the following sections. In session 2 we introduce 
general theoretical formalisms which we used in this paper. In section 3 we show and 
discuss our results obtained using the formalisms. In section 4 we summarize our results 
with some useful comments. 

2. Theoretical formalisms 

Theoretical analysis methods of measured optical data of s-wave and d-wave 
superconductors have been developed considerably [2, 4], One of important progresses, 
which has been successful for conventional superconductors [45, 46], is the boson- 
exchange model approach to study unconventional superconductors [2], One can extract 
the electron-boson spectral density function from measured optical self-energy [17, 18, 

47, 48]. The optical self-energy [S op (a;)] can be defined through an extended Drude 
model[6]. In the extended Drude model the scattering is inelastic (i.e. frequency 
dependent) while in the simple Drude model the scattering is clastic (i.e. frequency 
independent). When the scattering rate is frequency dependent its Kramers-Kronig 
counterpart is not zero any more. So we need a complex optical quantity to describe 
the inelastic scattering. The complex optical quantity is called the optical self-energy 
[E op (a;) = S° p (o;) -\-iT^{uj)\. The optical self-energy is a corresponding optical quantity 
of the quasiparticle self-energy [E 9P (a;)][9]. The imaginary part of the optical self-energy 
is related to the optical scattering rate [l/r op (a;) or the frequency dependent scattering 
rate]. The real part is related to the mass enhancement factor [A op (u;)] which is caused 
by the correlation among electrons [6]. 

Here we introduce formalisms which we used in this paper. The formalisms are a 
series of equations of the reverse process of the usual optical data analysis. We start 
from a model electron-boson spectral density function [/ 2 y(cu)] and use the Allen’s 
fromula[10, 15] to obtain the imaginary part of the optical self-energy The 

Allen’s formula is an integral equation which relates linearly the electron-boson spectral 
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density function [J 2 x(o;)] to the imaginary part of the optical self-energy [E^o;)]. The 
Allen’s formula can be described as follows: 


- 2E°» = 


T op (uj) 


dn r 2 x(n) k(u, o) + 


rZ P (u) 


( 1 ) 


where 1/t° p (o;) is the optical scattering rate, K(uj, hi) is the kernel for the Allen’s integral 
equation and 1 / T Zp( u ) * s the hnpurity scattering rate. We have different kernels for the 
normal state (here T — 0 case) and s- and d -wave superconducting states. The kernels 
can be written as follows: 
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where (~)(x) is the Heaviside step function (i.e. 1 for x > 0 and 0 for x < 0), E(x) is the 
complete elliptic integral of the second kind, and ■ -)g stands for the angular average 
over [0, 7 r/ 4 ]. A 0 and A 0 ( 6 *)[= A 0 008 ( 26 *)] are the superconducting gaps for s- and d- 
wave symmetry cases, respectively. 

The impurity scattering rates [1/r^ (a?)] for normal and s- and d- wave 
superconducting states can be described as follows[10, 15]: 
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where l/rf mp is the impurity scattering rate, which is a constant. We note that for the 
superconducting cases the impurity scattering rates are frequency dependent. 

In the extended Drude model formalism[ 6 , 36], which can describe correlated charge 
carriers, we can relate the optical conductivity, a(uj)[= oq(o;) + 202 ( 0 ;)] to the optical 
self-energy, E op (o;)[= E° p (o;) + iTff {co)\ as follows [ 6 ]: 
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where I2 p is the plasma frequency of charge carriers. The imaginary and real parts 
of the optical self-energy can be related to respectively the optical scattering rate 
[1 /t° p (o;)] and the optical mass enhancement factor [A op (a;)] as — VYffipj) = l/r op (u;) 
and — 2Tff{co) = \° p {co)oo = [m*(co)/m e — ljcu, where m*{co) and m e are the enhanced 
mass from the correlation and the electron mass, respectively. We note that the real 
and imaginary parts of the optical self-energy form a Kramers-Kronig pair. Therefore if 
you know one of them we can get the other using the Kramers-Kronig relation between 
them. 

In normal state at T = 0, —2X° p (a;) can be written analytically as follows [8, 10, 49]: 


- 2Y>° p (oo) = 2 / dm 2 X (n) 
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We also introduce an interesting relation at zero frequency as follows: 
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The second step in the first equation is obtained by applying the L’Hopitals’s rule. The 
last equation is called the coupling constant which we denote as A (see in Fig. 5) and 
oo c is the cutoff frequency. 

The dynamic dielectric function [e(u;) = ei(oo) -Me 2 (o;)] can be related to the optical 
conductivity as follows: 


47T _ 

e(co) = e H + i — a(oo) 

CO 


(7) 


where e# is the high frequency background dielectric constant. 

Finally one can calculate the single bounce reflectance [-R(cu)] at normal incidence 
from the dynamic dielectric function using the Fresnel equations and the relations 
between optical constants[5]. The reflectance can be described as follows: 


R(u) = 


N(co) - 1 
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where N is the complex index of refraction. N = n + in = \/f, where n and k are 
called the index of refraction and the extinction coefficient, respectively. We note that 
if a sample is not transparent its reflectance is the single bounce one. Here 1 is the 
dielectric constant of the vacuum or air. In the following section we will show and 
discuss about the results of our calculations using the formulas of the reverse process. 
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Figure 1 . (Color online) Two input model electron-boson spectral density functions 
(MMP and MMP+peak) in the upper frame, the calculated optical scattering rates 
for the three cases (one normal and two SC cases) in the lower frame, and impurity 
scattering rates for the three cases in the inset (see in the text). 

3. Results and discussions 

3.1. Model electron-boson spectral density, / 2 x(cu) 

For our calculations of the reverse process we started with two input model electron- 
boson spectral density functions, J 2 y;(u;), which are based on experimental results of 
hole doped cuprates[2]. One model function consists of a broad mode, of which shape 
is a phenomenological antiferromagnetic spin fluctuation proposed by Millis, Monicn, 
and Pines[35]. This is called the MMP model and described as / 2 x(tu) = 

s f 

where A s (= 50 meV) and u: s f (= 60 meV) are the amplitude and the characteristic 
frequency of the spin fluctuation mode, respectively. The MMP model has a maximum 
at uj s f and is displayed in the upper frame (red solid line) of Fig. 1. The other 
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model spectral density function consists of two components: the MMP and a relatively 
sharp Gaussian peak. The Gaussian peak is modeled as the magnetic resonance mode 
observed by an inelastic neutron scattering [24, 30] and is also called as the optical 
resonance mode[6, 18]. We denote this model as an MMP+peak model and describe as 


I\(u:) = 


A s Ld 

“ 2 +“i f 


+ 


V2ir{d/2. 35) 


exp 


(ui-Up ) 2 


, where A p (= 50 meV), c o p (= 30 meV), 


2(d/2.35) 2 

and d (= 10 meV) are the amplitude, the peak frequency, and the width of the Gaussian 
mode, respectively. The MMP+peak model is also displayed in the upper frame (blue 
dashed line) of Fig. 1 along with the MMP model (red solid line). 

The physical quantities related to the electron-boson spectral density function are 


the coupling constant [A = 2 f“ c I 2 x(k2)/QdQ, see also in Eq. (6)], the average electron- 
boson frequency ( u>i n ), and the maximum superconducting transition temperature 
(T '™ ax f The maximum superconducting transition temperature can be calculated in 
a generalized McMillan formalism [4, 23], which can be described as follows: 


k R Tf iax ^ 1.13 hu ln exp 


1 + Ai 
A . 


( 9 ) 


where ks is the Boltzmann constant, h is the reduced Planck’s constant, and uji n 
is the logarithmically averaged electron-boson frequency which is defined as c<y n = 
exp [(2/A) f Q ' In Q / 2 y (G)/G dtt\ where uj c is the cutoff frequency, we use u> c = 600 meV. 
In Table 1 we show the calculated values of the coupling constant, the average electron- 
boson frequency, and the maximum superconducting transition temperature. We note 
that the superconducting transition temperature is influenced mostly by the average 
frequency. We also note that these values in the table are comparable to reported ones 
for cuprates[17, 18]. 


Physical quantities 

A 

uiin (meV) 

T™ ax (K) 

MMP model 

2.424 

50.67 

161.8 

MMP+peak model 

5.830 

36.61 

148.8 


Table 1. The calculated coupling constant (A), the average frequency (wz„), and the 
maximum superconducting temperature (T™ ax ) of our two model input electron-boson 
spectral density function (see in the text). 


3.2. Reverse process 

Now we performed the reverse process starting from the two input electron-boson 
spectral density functions, I 2 x(oj). We obtained the imaginary part of the optical self¬ 
energy (or the optical scattering rate) for three cases using Eq. (1) and (2) and the 
input electron-boson spectral density functions. The three cases are a normal (N) state 
at T = 0 K and two superconducting (SC) states with s- and d- wave gap symmetries. 
The gaps at T = 0 K are A 0 (= 30 meV) and A 0 cos (2 6) for s- and d- wave SC cases, 
respectively. The obtained optical scattering rates [1/t° p (o;)] for the three cases using 
the two different I 2 x{ov) are displayed in the lower frame of Fig. 1. While the MMP 


















Reverse process of usual optical analysis of bo son-exchange superconductors: impurity effects 



E 

A 

j; 

cT 


"a, 

j; 

CN 

b 


Figure 2. (Color online) The real (lower left frame) and imaginary (upper left frame) 
parts of the optical self-energy for all three cases are shown in the left frames. The 
real parts are obtained from corresponding imaginary parts using a Kramers-Kronig 
relation. The real (upper right frame) and imaginary (lower right frame) parts of the 
optical conductivity are obtained using the extended Drude formalism for all three 
cases and displayed in the right frames. 


model gives gradual increases for all three cases the MMP+peak shows sharp rises with 
onset frequencies near the Gaussian peak energy (uj p = 30 meV) for the normal and the 
peak energy plus two times of the SC gap (u> p + 2A 0 = 90 meV) for the s- wave SC. For 
the d- wave SC case the sharp rise starts from near u p (30 meV) and ends at u p + 2A 0 
(i.e. 90 meV). We also calculated the impurity scattering rates [1/r,^ (a;)] for the three 
cases using Eq. (3) with l/rj mp = 1.0 meV and A 0 = 30 meV and displayed them in the 
inset of Fig. 1. While the normal impurity scattering rate is independent of frequency 
both SC impurity scattering rates show strong frequency dependencies below and above 
2 Ao and then become the normal impurity scattering rate value at high frequency region 
above 500 meV. 

Now we obtained the real part from the imaginary of the optical self-energy using 
a Kramers-Kronig relation. In Fig. 2 we display the obtained real parts of the optical 
self-energy [i.e. —2E° p (u;)] in the lower left frame for the three cases with two input 
model I 2 x{u) along with their corresponding imaginary parts [i.e. — 2 T, 2 P (co)} in the 
upper left frame. The results of normal state agree with the reported ones which were 
obtained in a similar situation [50]. For two SC states with the MMP+peak / 2 y(o;) 
there are prominent peaks near uj p + 2A 0 (i.e. 90 meV) where the sharp rises are in 
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Figure 3. (Color online) The real and imaginary parts of the optical self-energy with 
two representative impurity scattering rates (l/ri mp ), 1.0 meV (a clean limit) and 30.0 
meV (a dirty limit). In the insets magnified views of the optical self-energy for the 
dirty limit case in low energy region are displayed. The red arrow indicates a sharp 
peak in the real part of the optical self-energy which is induced by the impurities. 


the optical scattering rates. These peaks are originated from the sharp Gaussian peak 
in the input / 2 x(w). In the corresponding real part of the normal state we observe a 
kink-like feature (a abrupt slope change) near oj p [i.e. 30 meV); the curve would have a 
logarithmic singularity in the slope at c o p if we were using a pure delta function instead 
of the sharp Gaussian function[8] with a finite width. The characteristic features of 
I 2 x(p}) are clearly be seen in the optical self-energy. 

Furthermore in the extended Drude formalism [see Eq. (5)] we obtained the optical 
conductivity with the plasma frequency (D p ), 2.0 eV (or = 16,130 cm -1 ) and displayed 
the optical conductivity spectra for the three cases in the right frames [<7i(a;) in the 
upper and a 2 (cj) in the lower] of Fig. 2. For this impurity-free [i.e. l/T irnp = 0) case we 
can see only the Holstein boson-assisted absorption piece in the real part of the optical 
conductivity, which is a fraction [A/(l + A)] of the total spectral weight, where A is the 
coupling constant (see Table 1). The total spectral weight [W Sjto tai) is 7 t/120 times the 
plasma frequency squared, i.e. W S)to tai = f^+ Gi(oj) du = (7r/120) D 2 , where G\ (uj) is in 
units of n -1 cm -1 and Q p in units of cm -1 . The total spectral weights for all three cases 
should be the same since we used the same plasma frequency for the three cases in the 
extended Drude model. We compared the incoherent Holstein boson-assisted absorption 
spectral weights (= 4.6 x 10 6 and 5.6 x 10 6 !7 -1 cm -2 obtained from the integration of the 
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Figure 4. (Color online) Calculated real parts of the optical conductivity [cri(w)] and 
real parts of the dynamic dielectric function [ei(w)] times -oj 2 for all three cases with 
two different input / 2 y(w) with two representative impurity scattering rates (l/ri mp ), 
1.0 meV and 30.0 meV. We also show the real part of the optical conductivity spectra 
of the Drude modes (black dash-dotted lines) with corresponding impurity scattering 
rates. 

optical conductivity for the MMP and the MMP+peak, respectively) with the fractions 
[A/(l +A) = 0.708 and 0.854 for the MMP and the MMP+peak, respectively] of the total 
spectral weight (6.7x 10 6 0 -1 cm -2 ), = 4.7x 10 6 and 5.7x 10 6 f2 _1 cm~ 2 for the MMP and 
the MMP+peak, respectively, ffere we took the coupling constants (A) for the MMP 
and MMP+peak from Table 1. The resulting incoherent spectral weights obtained from 
two different ways show good agreements for the both input I 2 x(u) cases. In cr 2 (u;) 
we can see that the sharp Gaussian peak in the MMP+peak model appears as very 
prominent kinks (abrupt slope changes) near uj p + 2A 0 (i.e. 90 meV) for two SC states 
and a kink near oj p (i.e. 30 meV) for the normal state. At the same frequencies where 
the kinks are we can see sharp rises in the real part of the optical conductivity, ai(to). So 
we can clearly observe the characteristic features of / 2 y ( oj) in the optical conductivity. 

3.3. Impurity effects 

Now we added various impurity scattering rates (1, 5, 15, and 30 meV) to the impurity- 
free optical scattering rates as in Eq. (1), (2) and (3) and used the same Kramers- 
Kronig relation to obtain the real parts of the optical self-energy. In Fig. 3 we display 
the optical self-energy for two representative impurity scattering rates, l/rj mp = 1.0 
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Figure 5. (Color online) The partial sum [VF s (u;)] and the difference between normal 
and superconducting states, W Si n(u}) — W St sc{u), for all three cases with two different 
input / 2 x(w) with three representative impurity scattering rates (l/r,: mp ). 0.0, 1.0 and 
30.0 meV. We also show the partial sums of the Drude modes (black dash-dotted lines) 
with corresponding impurity scattering rate. We note that for 1 /Ti mp = 0 meV there 
is no spectral weight at finite frequency region. 


and 30.0 meV, which can be the clean (l/rj mp -C 2A 0 ) and the dirty (1 /r imp < 2A 0 ) 
limit for s-wave and approximately d- wave superconductors, respectively. We observe 
clearly a very sharp peak (marked with a red arrow in the inset of the right lower 
frame) in the real part of the optical self-energy for the dirty limit case (1 /Ti mp = 
30.0 meV) at 2A 0 (he. 60 meV) for s-wave SC case. At the same energy there is a 
sharp rise in the imaginary part of the optical self-energy as shown in the inset of the 
right upper frame. The sharp features seem to exist in the d- wave SC case but they 
are significantly suppressed; we can see it if we look at very carefully near 2Ao (he. 60 
meV). This sharp feature is getting more pronounced as the impurity level increases, he. 
this feature can be induced by impurities. The sharp feature may not be observed in d- 
wave superconductors by infrared spectroscopic technique, which measures the averaged 
response over the anisotropic Fermi surface (or the k space). But a similar impurity- 
induced sharp feature was observed in B^S^CaC^Og+s (a d- wave superconductor) 
by high-resolution ARPES[28, 51, 52]. This feature showed strong temperature and 
momentum dependencies. It also was interpreted theoretically as a result due to forward 
scattering in the d- wave cuprates [53]. 

We also calculate the optical conductivity from the optical self-energy using the 
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extended Drude formula [see Eq. (5)] with the same plasma frequency (it p = 2.0 eV) 
at various impurity scattering rates (1, 5, 15, and 30 meV). In Fig. 4 we display 
the calculated real parts of the optical conductivity and the real parts of the dynamic 
dielectric function for two representative impurity scattering rates, 1/Tj mp = 1.0 and 
30.0 meV. We also display the real parts (dot-dashed black lines) of the corresponding 
Drude optical conductivity in the upper frames. The Drude conductivity is completely 
coherent; there is no incoherent Holstein boson-assisted absorption since the charge 
carriers are essential free (or A = 0). In other words the optical self-energy is zero, 
i.e. E op (a;) = 0, only elastic impurity scattering exists. For the normal cases (red solid 
and red dashed lines) we can see that the additional spectral weights appear at low 
frequency region in the optical conductivity compared with corresponding spectra for 
1 /^imp = 0 (see in the upper right frame of Fig. 2). Now the coherent portion of the 
total spectral weight [i.e. W S:to tai/( 1 + A)] also shows up in the conductivity spectra in 
a finite frequency range; the coherent spectral weight was confined at zero frequency for 
1 /^imp = 0. For 1/Tjmp = 30.0 meV (or the dirty limit) case we can see the gap features 
(i.e. sharp rises) clearly at 2A 0 (i.e. 60 meV) in ai(u>) of the s-wave superconductors. 

For the d- wave SC cases we also can see the gap features (i.e. sharp rises) at 2A 0 but 
the features are suppressed significantly by the anisotropic SC gap. For l/rj mp = 1.0 
meV (or the clean limit) case those SC gap features cannot be seen very clearly. In the 
optical conductivity calculated by the reverse process we can observe the Gaussian peak 
in I‘ 2 x(uj) which appears as a dip near uj p + 2A 0 (i.e. 90 meV) for both SC cases and 
near uj p (i.e. 30 meV) for the normal case. 

In the lower frame of Fig 4 we display — ei(uj)uo 2 for all three cases. In principle 
this quantity at zero frequency gives the superfluid plasma frequency squared (D^ p ), i.e. 
lirn^of—ei(o;)a; 2 ] = f^ p . We clearly can see the superfluid plasma frequency decreases 
as the impurity scattering rate increases for both superconducting cases; this result 
is expected for a superconductor moving into the dirty limit [37, 38]. The resulting 
superfluid plasma frequencies for various levels of impurities are calculated and shown 
in Fig. 6. For the normal case this quantity should be zero i.e. lim^of - dO^V 2 ] = 0 
since there is no superfluid current. In —ei(co)co 2 the SC gap feature appears as a very 
sharp dip at 2A 0 (i.e. 60 meV). We also can clearly see the sharp Gaussian peak 
feature in the electron-boson spectral function near uj p + 2Ao (i.e. 90 meV) for both 
superconducting states and near u p (i.e. 30 meV) for the normal state, which appear 
as broad valleys. 

Now we extract the superfluid plasma frequency with another method. In this 
method we monitor the spectral weight redistribution (or the missing spectral weight) 
when the system experiences a phase transition from normal to superconducting states. 
This method was originated from the Ferrel-Glover-Tinkham (FGT) sum rulc[40, 41], 

We denote the superfluid plasma frequency extracted using this method as Q sp ,fgt and 
call as the FGT superfluid plasma frequency. We can describe fl sp ,FGT in the FGT sum 


Reverse process of usual optical analysis of bo son-exchange superconductors: impurity effects 



0 5 10 15 20 25 30 

1/t. (meV) 

imp v ' 


Figure 6. (Color online) Impurity dependent superfluid plasma frequencies obtained 
by two different methods (see in the text) for two different superconducting gap types 
(s-wave: upper frame and d- wave: lower frame) with the MMP and MMP+peak input 
i 2 xM- 


rule formalism as follows: 

120 

Qs P ,fgt(u) 2 =-[Ws.jvO^O — Ws,sc( w )] (10) 

7T 

where W s ^n(lo) and W S) sc{u) are the spectral weights of normal and superconducting 
states, respectively. We note that is in units of cm -1 and W s (uj) in units of 

n -1 cm -2 . For both normal and superconducting states the spectral weight is defined as 
W s (lo) = J Q + cxi(a/) duj'. In the upper frames of Fig. 5 we display the spectral weights 
( W s {uj )) for one normal and the two SC cases with two input model / 2 x(u;) and three 
representative impurity scattering rates (0.0, 1.0, and 30.0 meV) in an extended spectral 
range (up to 6.0 eV). The corresponding ^tfp^GT^) f° r the three representative impurity 
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Figure 7. (Color online) Calculated reflectance spectra for all three cases with two 
different input / 2 x(w) with two representative impurity scattering rates (1/Ti mp ), 1.0 
and 30.0 meV. We also show reflectance spectra of the Drude modes (black dash-dotted 
lines) with corresponding impurity scattering rates. In the insets we show magnified 
views of the reflectance in low frequency region. 


scattering rates are calculated and displayed in the lower frames. Interestingly, we have 
a finite superfluid plasmas frequency for 1/T imp = 0.0 meV case, which indicates that 
some portion of the incoherent Holstein boson-assisted absorption is condensed into the 
superfluid state. At the same time we also can see that the total coherent portion of 
the spectral weight is smaller than the superfluid spectral weight for l/rj mp = 0.0 meV 
case, i.e. W s>total [1/(1 + A)] < W SyFGT = (7 t/120 ) ^% yFGT for both input J 2 y( u). This 
supports the idea that some incoherent Holstein spectral weight was condensed into the 
superfluid density. Since there is an overshot in the FGT superfluid plasma frequency 
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Figure 8. (Color online) Comparison of the optical scattering rate (1 /t° p (u>)) and 
1 - R(u). 


which extends up to around 2000 cm -1 a reasonable estimate for the cutoff frequency 
for the FGT plasma frequency can be around 2000 cm -1 . 

In Fig. 6 we display the superfluid plasma frequencies obtained with two different 
methods as a function of the impurity scattering rate. The plasma frequencies obtained 
using Vt 2 p = lim^of—< e i( UJ ) ^ 2 ] are shown with solid symbols and those obtained using 
^sp,fgt{ u ) = ^p\W s ,n{u) — Ws’Sci 1 ^)] are displayed with open symbols. In low impurity 
scattering rates (or in the clean limit) below l/rj mp = 5.0 meV n^ FGT (w) decreases very 
quickly with reducing the impurity scattering rate because the coherent portion is too 
narrow to be taken into account completely. This means that to get a good agreement 
between the two methods the system should contain enough level of impurities in order 
to make the coherent portion broad enough. But in high impurity scattering rates (or in 
the dirty limit) above 15.0 meV results from both methods agree each other quite well. 
We observe that the superfluid plasma frequency depends on the impurity scattering 
rate; as the impurity scattering rate increases the superfluid plasma frequency gradually 
decreases. The difference between the two superfluid plasma frequencies obtained by 
the two methods can be considered as an uncertainty for the FGT superfluid plasma 
frequency. We also observe that the superfluid plasma frequency of a system with 
the MMP+peak in I 2 x{u) shows more robust against the impurities than that with 
the MMP alone in I 2 x{uj). Interestingly, the superfluid plasma frequencies of s-wave 
superconductors seem to be more robust against the impurities than d -wave ones. 

Eventually, we calculated reflectance spectra [-R(w)] using the Fresnel equation 
with the dynamic dielectric functions [see Eq. (8)] with €h = 1-0. We note that 
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the calculated reflectance spectrum is at normal incidence. In Fig. 7 we display the 
calculated reflectance spectra for two representative impurity scattering rates (1 /r imp = 

1.0 and 30.0 meV). We also display Drude reflectance spectra (dash-dotted lines) with 
corresponding impurity scattering rates. We can see that correlated electron systems 
show much suppressed reflectance below the plasma frequency (h2 p = 2.0 eV) compared 
with the free electron (or simple Drude) system. In the insets we show magnified views 
of reflectance spectra in low frequency region. For l/rj mp = 30.0 meV the reflectance 
spectra of the s-wave SC case show clearly the superconducting gap feature at 2A 0 
(i.e. 60 meV); below the energy the reflectance becomes essentially 1.0. For the d- wave 
SC case the gap feature is significantly suppressed; the reflectance becomes 1.0 at oj = 

0. We also can see characteristic features which are caused by the input / 2 y(o;). To 
see those features more clearly we compare reflectance spectra with the corresponding 
optical scattering rates; we display 1 — R(uj) and l/r op (u>) in Fig. 8. Those two spectra 
look quite similar to each other. To show the similarity explicitly is not trivial since the 
reflectance is a quite complicated quantity theoretically. However the similarity can be 
understood roughly, as follows: in general at low frequency region a low scattering rate 
gives high conductivity which results in high R(oj) (or low 1 — R{w)) and vice versa for 
a high scattering rate. So roughly one can say that 1/t° p (lo) oc 1 — R(oj). Therefore one 
may be able to see characteristic features of J 2 y(o;) as well as the gap features in the 
raw measured reflectance spectrum as one can see them in the optical scattering rate. 

We note that thermal excitations may smear the distinction of the features and/or shift 
the position of the SC gap features. 

3-4- Comparison with experimental data 

Now we compared measured experimental data with the corresponding optical data 
obtained by the reverse process. In Fig. 9 we display measured reflectance data of a 
d- wave superconductor, optimally doped Bi 2 Sr 2 CaCu 20 8+( 5 (Bi-2212) with T c = 96 K, 
at T = 27 K (superconducting state) and 102 K (normal state) and optical constants 
obtained using the usual optical data analysis [2 7] in (a-f) frames and the corresponding 
optical constants including reflectance obtained from the reverse process for d- wave SC 
case with the MMP+peak model / 2 y(cu) for l/r; mp = 1.0 and 30.0 meV in (g-1) frames. 

In (b) and (h) frames we can see that the results with 1/T irnp = 30.0 meV show better 
agreements to the experimental date, which indicates that certain impurities caused by 
defects and/or thermal excitations exist in the experimental data. The experimental 
data would show temperature smearing effects from thermal excitations. Overall all 
optical data displayed show quite good agreements even though experimental data are 
at finite temperatures while our calculations are done at zero temperature. We note that 
the electron-boson spectral function of Bi-2212 shows an additional intrinsic temperature 
dependent evolution as well [2, 18]. We also note that in the calculations we used the 
same input / 2 y (o;) for both normal and superconducting states and for l/r; mp = 1.0 and 
30.0 meV. We note that the real part of the optical self-energy shows relatively small 
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Figure 9. (Color online) Comparison of (a-f) measured reflectance and other optical 
constants obtained using the usual analysis of an optimally doped B^S^CaC^Og+a 
(Bi-2212) and (g-1) results obtained from the reverse process for d -wave SC case with 
the MMP+peak model / 2 \(u;) for l/T imp = 1.0 (solid lines) and 30.0 meV (dotted 
lines). 


impurity dependence. Actually for the normal case — 2Y,° p (uj) should not depend on the 
impurity level at all since a frequency independent impurity scattering rate is added and 
its Kramers-Kronig counterpart is essentially zero. However for superconducting cases 
since the impurity scattering rates show frequency dependence near 2Ao their Kramers- 
Kronig counterparts are not zero any more. This gives the difference between —2S° p (tu) 
obtained with different impurity scattering rates near 2Ao- While the positions of the 
peaks of calculated and measured — 2Y,° p (u}) have similar energies (he. 100 meV in 
frame (e) and 90 meV in frame (k)) the peaks of two (experimentally obtained and the 
input) / 2 y(cu) are located at quite different energies (he. 60 meV in frame (f) and 30 
meV in frame (1)). This difference comes from different superconducting energy gaps 
(he. Ao — 20 meV for experimental data and Ao = 30 meV for the calculated ones). 
As we mentioned in section 3.2 the peak in —2Tff{uj) is located at uj p + 2A 0 . So we can 
understand that 100 meV ~ 60 meV + 2x20 meV for the peak in experimental data 
(frame (e)) and 90 meV = 30 meV + 2x30 meV for the peak in calculated one (frame 
(k)). These good agreements confirm that our reverse process is quite reliable. 
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4. Conclusions 


In this study we have performed a reverse process of the usual process for optical 
data analysis. The reverse process is as follows: We started with the input model 
electron-boson spectral function, / 2 x(u;). We obtained the imaginary part of the 
optical self-energy using Allen’s formulas, the real part of the optical self-energy 
using a Kramers-Kronig relation, the optical conductivity using the extended Drude 
model, other optical constants using well-known relations between optical constants, 
and eventually reflectance using the Fresnel equations. We applied this reverse process 
to three cases (one normal and s- and d- superconducting cases with the superconducting 
gap Ao = 30 meV) with two (MMP and MMP+peak) input model / 2 y(u;) which 
are based on the reported experimental results. Since impurities are not avoidable 
experimentally we included various levels of impurities (from the clean to the dirty limit: 
1 / T imp = 0 ~ 30.0 meV) in the optical scattering rate and performed the same reverse 
process to investigate impurity dependent optical properties. In the clean limit (for 
example, l/ri mp = 1.0 meV) for the both superconducting cases we do not observe the 
superconducting gaps clearly but still the characteristic feature (i.e. the sharp Gaussian 
peak and the MMP model) of / 2 x(o;) appears definitely in all optical constants including 
reflectance. For the normal case we also are able to see the feature of J 2 x(o;) in the 
optical constants including reflectance. In the dirty limit (for example, l/rj mp = 30.0 
meV) the superconducting gap appears definitely in s-wave superconductors but the 
gap feature is suppressed significantly in the d- wave one. In d- wave superconductors 
the optical constants are averaged response signals over the anisotropic Fermi surface 
which causes the suppression of those features. Furthermore in the dirty limit the 
characteristic features (the Gaussian peak and the MMP model) of J 2 x(o;) also can be 
seen clearly in the optical constants and measured reflectance spectrum. We also find 
that in the clean limit the superfluid density obtained by the FGT sum rule cannot 
capture whole coherent electrons which participate in the superconductivity because 
the coherent absorption peak is too narrow. To get accurate superfluid density the 
SC system should be in a dirty limit. Our study shows that the superfluid density 
obtained by the FGT sum rule becomes quite accurate when 1/T imp is large enough, 
above 15 meV, which is comparable to the superconducting gap, A 0 (30.0 meV). We 
observe that the superfluid density decreases as the impurity level increases for both 
SC cases. We also observe that the s-wave superconductor is more robust against the 
impurities than the d- wave one. Our results will help to understand the usual optical 
analysis process thoroughly as well. Since the optical spectra show clearly the mediated 
boson features of the boson-exchange superconductors collecting experimental data at 
various (temperature, doping, magnetic field and so on) conditions and obtaining J 2 y( u) 
applying the usual process can be a first step to reveal the microscopic origin of the 
exchange boson. 
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